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1. The project, in a nutshell 

Two coupling constants, F and £, parameterize the leading-order chiral effective Lagrangian 

J^eff = —Tr(d^Ud^) - -TrL#(J/ + [/+)]. (1.1) 

A long-standing problem for lattice simulations is to determine them directly from QCD. Typically, 
this is done by fitting a correlation function to some theoretical formula parameterized by F, £ and 
m % . Usually one has to measure a correlation function of operators at long distance to be sensitive 
to these parameters. This makes the measurement expensive. In addition, the simulations do 
not give F and £ directly; instead, they give the mass-dependent decay constant and condensate, 
which depend, through the standard formulas of chiral perturbation theory, on F and £ and on the 
coefficients of higher-order terms in Jz? '. So the determination of F and £ from the lattice is rather 
indirect. 

A series of papers by Akemann, Damgaard, Heller, Osborn, Splittorf, Svetitsky, and Toublan^], 
IH D, ||] provide an alternative path toward a direct measurement of F and £. This is done via the 
properties of eigenvalues of the Dirac operator, from simulations in the epsilon regime {m n L <C 1 
but AL S> 1 for box length L and A any nonperturbative parameter of QCD except m n ). The par- 
ticular implementation we used is that of Ref. (g). An imaginary isospin chemical potential /I is 
coupled to a doublet of quenched quarks. The correlator of the eigenvalues from an ordinary /j. = 
simulation, A,-, with the quenched eigenvalues, Ay, 

P ^ onn (x,y) = (Z8(x-X i )^8(y-Xj))-(^8(x-X i ))(^S(y-lj)) (1.2) 

' j ' j 

is a function of £ and F which is given by Random Matrix Theory (RMT). The particular formula 
is given by Eq. (3.49) of Ref. |^|, in terms of the usual rescaled variables X{LV , XjLV, the rescaled 
mass m q LV, and the rescaled isospin chemical potential 8 = pLF^/V. V is the volume. 
To combine data sets, we integrate the data. Refs. [Q] and [^] suggest using 

f Cmax 

C(x,Cmax)= dyp(x + y,y) . (1.3) 
Jo 

We did one more integral: 

I(X,Uax)= [_ C(x^ max )dx. (1.4) 

The first integral C(x) shows a spike at x = whose width goes roughly as S 2 . Then I(X) will show 
a sharp step at X = 0. 

In our numerical simulations we used overlap fermions|§, They preserve the full SU (Nf) (8) 
SU (Nf) chiral symmetry (including the anomalous singlet current and the index theorem). Our 
data set uses a lattice volume of 12 4 points. The lattice spacing a, as determined from the Sommer 
parameter ro [§], is ro/a = 3.71(5). The pseudoscalar mass is am n = 0.329(3). This is not in 
the epsilon regime, but experience shows RMT is robust enough to work well outside the epsilon 
regime, and we are making a pilot study. In our conversions to physical units we use ro = 0.5 fm. 



Details of our action are to be found in Refs. M, 10, 11]. The one new ingredient is the gauge 
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Figure 1: Data (squares) and fit from |v| = 1. A cut C, max = 0.07 is enforced. The best fit values for this 
bootstrap sample are LlV = 147 and aFi = 0.071. 
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Figure 2: Best-fit values of qFl, varying £„, 



connection: we used the differentiable hypercubic smeared link of Ref. Q12|]. We collected about 
400 thermalized HMC uajectories of unit length and analyzed 30 lattices in topological charge 
sector v = and 75 |v| = 1 ones. To do the fitting, we took bootstrap averages of the integrated 
data. 

We reproduce two figures from the paper: Fig. [l] shows an example of a fit from one of our 
bootstraps. Fig. ^| shows the dependence of fits on C, max - These fits all include only the three lowest 
(but nonzero) eigenvalues in a topological sector. Smaller £ mar means that less data is included in 
the fit; larger ^ max means that more eigenmodes are needed to saturate the correlation function. It 
appears that results for this data set do not depend on the choice of C, max . 

In finite volume, the eigenvalue correlator is not a direct function of £ and F, but of the 
quantities Ll and Fi which reflect the finite volume. For the condensate, this dependence is known: 
Ll = Pz£ where 



(1.5) 
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with A(0) the contribution to the tadpole graph (propagator at zero separation) from finite-volume 



image terms. In the epsilon regime, A(0) = —Dj\jV and D depends on the geometry Q13Q. (It is 
0.1405 for hypercubes.) For N f flavors, C z = -(Nf- l)/N f . 

A second complication for the condensate is that £ is scheme-dependent. We perform the 
conversion from lattice regularization the MS using a measurement of the matching factor from the 
RI-MOM scheme as an intermediary. Correcting for finite volume and converting schemes, we find 

r L(MS, jU = 2 GeV) 1/3 =0.594(13) (1.6) 

or 

(L(MS,H = 2 GeV)) 1 / 3 = 234(4)MeV (1.7) 

This is consistent with an analysis of cumulants of eigenvalues, and with other determinations of 
the N f = 2 condensate[0, [B|, 



F is also expected to be volume dependent, and the correction is also likely to be of the form 
Fl = PfF with 

Pf = 1+C F ^A(0) + ... (1.8) 
Unfortunately, there is no calculation (yet) of Cp. We found 

fo^L = 0.255(13) (1.9) 

or with ro = 0.5 fm, 

Fl = 101(6)MeV. (1.10) 

We expect that Cp is a number on the order 1-2, probably Nf dependent; that would be sufficient to 
lower F from Fl by 20 per cent or so and bring it in line with phenomenological estimates of about 



86 MeV [18, 1 



So at this point the calculation is incomplete: one needs the finite volume correction. However, 
the size of the error is very interesting: while typical large scale simulations can get f n to a fraction 
of a per cent, F is harder: for example, MILCpO]] quotes (their Table IV, Fit A, our conversion) 
r\F = 0.131(10) or F = 82(6) MeV, compared to their number for f K , which has a 0.3 per cent 
error ||2l|]). Our quite competitive error on F, with the expenditure of a tiny amount of computer 
resources, is what makes this calculation worth repeating. 

2. Questions from the conference 

• You used fat links. Aren't you worried they will mess something up? 

On the contrary. Everyone who has tried using fat links for light fermions has only good things 
to say about them, from the point of view of perturbation theory, exceptional configurations (for 
Wilson fermions), taste symmetry restoration (for staggered fermions) and efficiency in computing 
the overlap action. The one exception was MILC, who long ago[[22"|] tried to use APE smearing 
with 0.45 x 10 (!!!) smearing steps in a calculation of the properties of heavy-light mesons. This 
gave big scaling violations in fp,: too much smearing decoupled the quarks from the short distance 
part of the one gluon exchange potential. We are using hypercubic smearing. This is even more 
local than AsqTad - in fact, you cannot be more local. Readers might remember that the Alpha 



collaboration studied static-light properties with hypercubic links, with great success [23]. 
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• Why did you measure Fl Isn't f n more interesting? 



Actually, what is interesting is the ratio fic/fx, which can be used to determine the Cabibbo 



angle[|24|]. But then you have to ask, how do you measure f K l This is not done with simulations 
at the physical quark masses; instead, one performs simulations at any light quark masses (heavy 
enough to be cheaply simulated, light enough to be in the region of validity of chiral perturbation 
theory). Then the mass-dependence of parameters is fit to a chiral perturbation theory formula. 
The uncertainties between F, £ and the higher order terms cancel at the input quark mass values, 
for the particular observables which are being fit. However, for other observables (predictions) the 
larger error on F can be problematic. That's why this particular measurement scheme - provided 
the finite volume correction can be controlled - is so interesting. 

• Isn't there a problem using individual eigenvalues to determine F and £, rather than using 
the complete spectral density p(£)? 

We don't think so (and anyway, we used formulas which involved integrated eigenvalue distri- 



butions). Ref. [ )25[ ] showed that the spectral density and its moments, and even eigenvalues, renor- 
malize pretty much as one would naively expect. Notice that at low values of £, p(Q coincides 
with the distribution of the lowest eigenvalue pi(Q (because that is all that there is)! Akemann 
and Damgaard [^] have derived exact expression relating the individual eigenvalue distributions 
to integrals over the eigenvalue density and all the spectral correlators. If these are true, it is hard to 
imagine that something different would happen to the individual eigenvalue distributions and to the 
integrated ones. (They are also working on a single-eigenvalue version of the correlator we studied 



hereg7|].) 

• I don't have a chiral action. How much of this project could I do? 

Doing exactly what we did with fermions which do not recognize the index theorem could be 
tricky, since you have to work in sectors of fixed topology, and you have to project the nonzero 
eigenvalues, which are complex for a generic non-chiral Wilson-type fermion, onto the imaginary 
axis. This might be more trouble than it is worth. You might try looking at quantities in the 



epsilon regime which average over topology; the original papers by Hansen and Leutwyler[28] 
give formulas for the pseudoscalar and axial correlators, which can also be analyzed to get E and 
F. 

• Tell me more about the "Primme" package. 



The "Primme" package of McCombs and Stathopolous[g9|] is a library for computing eigen- 
values of Hermitian matrices. (For overlap fermions, that is all you need.) It implements a variety 
of algorithms, preconditioning and optimization options which can be used via a very convenient 



interface. It is about a factor of three faster than our version of the Bunk-Kalkreuter-Simma [ |30| ] 
Conjugate Gradient algorithm. 

• Can you suggest a good book about lattice QCD? 

Yes pi]. 
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